Frustrated classical Heisenberg model in 1 dimension with added nearest-neighbor 

biquadratic exchange interactions 
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The ground state phase diagram is determined for the frustrated classical Heisenberg chain 
with added nearest-neighbor biquadratic exchange interactions. There appear ferromagnetic, 
incommensurate-spiral, and up-up-down-down phases; a lock-in transition occurs at the spiral 
boundary. The model contains an isotropic version of the ANNNI model; it is also closely related 
to a model proposed for some manganites. The Luttinger-Tisza method is not obviously useful due 
to the non-linear weak-constraint problem; however the ground state is obtained analytically by the 
exact cluster method of Lyons and Kaplan. The results are compared to the model of Thorpe and 
Blume, where the Heisenberg part of the energy is not frustrated. 

PACS numbers: 75.10.Hk,75.30.Kz,75.47.Lx 



The ANNNI (antiferromagnetic next-nearest-neighbor 
Ising) model [J, HJ has Ising spins Sf, i.e. 2- valued ob- 
jects, located at points i on a simple cubic lattice with 
nearest-neighbor ferromagnetic interactions J± plus next- 
nearest-neighbor antiferromagnetic interactions Ji along 
one of the cubic directions, say x. Its ground state is 
the same as that of the linear chain (translated to all 
x-chains), whose Hamiltonian is 
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■h ^ s. 
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(i) 



n running over the integers, either — oo to oo or with peri- 
odic boundary conditions; also J\ < 0, J2 > 0, the latter 
embodying frustration or competition between the two 
terms. The ground state phase diagram, which depends 
only on 7 = J2/IJ1I, has a very simple structure: there 
is ferromagnetic orderingfor 7 < 1/2, up-up-down-down 
ordering for 7 > 1/2. 

The common isotropic version of ((T|) is the Heisenberg 
model obtained from dXJ) by the replacement S^S^ — > 
S ra • S m , in the classical version of which the spins are 
classical unit vectors. We will consider the classical ver- 
sion, which is in fact the mean field approximation to the 
quantum model Q. The corresponding ground state is 
ferromagnetic for 7 < 1/4, spiral for 7 > 1/4, the wave 
vector q varying continuously from as 7 increases past 
1/4. 

Thus, not surprisingly, there is great qualitative dif- 
ference between the (anisotropic) Ising and (isotropic) 
Heisenberg cases. A common way of interpolating be- 
tween these models is to consider the "XXZ" Hamilto- 
nian, obtained from the Heisenberg case by the replace- 
ment Jnm^n ' S m ► J nm (S n S m + S^S^) + J nm S n S m . 

This is anisotropic unless Jf lm = J r z lm in general. I want 
to consider a different connection, which maintains full 
isotropy, but nevertheless retains some characteristics of 
the Ising case. Namely, add a biquadratic term to the 



Heisenberg model: 

H = ^^[JiS„-S„ + i+ J2S„-S n+ 2— a(S„-S„ + i) 2 ], Sf = 1. 

(2) 

This, with J\ < 0, J2 > as above, is the model that will 
be addressed subsequently. That one may expect Ising- 
related ordering for large positive a can be anticipated 
because for J\ = J 2 = 0, the set of ground states is the 
set of collinear states, although there is degeneracy as to 
which axis all spins are parallel. Thus the entropy per 
spin is In 2 in the thermodynamic limit (the contribution 
of this rotational degeneracy disappears in the T.L.), the 
same as for the (non-interacting) Ising model. 

Biquadratic exchange has a long history of being found 
to be important in certain circumstances. E.g., one of the 
earliest works indicating appreciable effect of such inter- 
actions is in the paramagnetic resonance experiments of 
Harris and Owen [5(, that studied the nearest-neighbor- 
pair spectrum of Mn 2+ ions in MgO. They find that a 
value of j = 0.05Jin the Hamiltonian JS a -S{,— j(S a -S{,) 2 
gives a much improved and rather good fit to their mea- 
surements. The assumption that the coefficient 0.05 indi- 
cates a small effect would be wrong: In fact the correction 
to the Heisenberg term is almost a factor of 2 (i.e. 100%) 
for some of the observed and calculated Lande intervals; 
this comes from the large spin factors involved. The mi- 
croscopic origin and an order-of-magnitude estimate were 
discussed by Anderson. Q For more recent work see [?J 
and references therein, and below. I note, in particular, 
the consideration by Thorpe and Blume [8| of the special 
case of ©, J2 = 0. 

The well-known Luttinger-Tisza method appears to be 
not useful for finding the ground state of ([2]) because of 
the non-linearity introduced into the equations for sta- 
tionarity of H subject to the weak constraint, 
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Instead I turn to the rather unknown cluster method of 
Lyons and Kaplan which is tractable and solves the 
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FIG. 1: Angles made by the spins in a cluster of 3. 



problem rigorously. Briefly recall that method. Assume 
periodic boundary conditions. Then one easily verifies 
that @ can be rewritten as 



H 



E 



h c (Si-i, Si, Sj-|-i) 



(3) 



where the "cluster energy" 
/i c (Si,S 2 ,S 3 ) = 

i{Ji(Sr .S 2 + S 2 .S 3 ) 
- a [(S 1 .S 2 ) 2 + (S 2 -S 3 ) 2 ]} + J 2 S 1 -S 3 
involves 3 neighboring spins. Clearly 

H > V"min/i c (Si_i,S 4 ,S i+ i). 



(4) 



(5) 



One can easily find the minimum of h c . If the correspond- 
ing state "propagates" , i.e. if there is a state of the whole 
system such that every set of 3 successive spins gives the 
minimum h c , then according to ([5]), this state will be a 
ground state of H. This is the LK cluster method as ap- 
plied to the present problem. The method is not limited 
to 1 dimension or to periodic Hamiltonians. [|| 

Now let's minimize h c . First consider coplanar states, 
and label the angles made by the end spins with the 
central spin 8,6', assumed with no loss of generality to 
be up, as shown in Fig. 1. The cluster energy is then 

h c (9,9') = -i(cos0+cos6>')+7Cos(6>-6>')--(cos 2 6»+cos 2 ( 

2 (6) 

where for simplicity I have put J\ = — 1 (ferromagnetic), 
and used the previous definition 7 = J 2 /| J\\. Differenti- 
ating gives the conditions for stationarity 
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Solutions are 



— 7sm( 
7 sin(6* 



9') + asin9cos9 



(7) 



(9,9') 
(9,9') 

cos 6*o 



(0, 0), (0, tt), (tt, 0), (tt, tt) (Ising type), and 
(9q, —9q) (spiral type), where 
1 



2(2 7 - a) 



for |2(27-o)| > 1. 



(8) 



The (tt, tt) solution (which leads to the ordinary anti- 
ferromagnetic state) is never lowest because we have as- 
sumed Ji < 0. The (0,0) solution obviously propagates 
as the ferromagnetic state. The solutions (71", 0), (0, tt), 
i.e. (|, I, f), (1, 1, |) plus their degenerate reversed spin 
counterparts can easily be seen to propagate in the up- 
up-down-down state. The solution (9q, — #o)> degener- 
ate with its uniform rotations, obviously propagates in a 
simple spiral 



S„ = x cos nQo + y sin u9q , 



(9) 



x,y being any pair of orthonormal vectors. Such states 
were first discussed long ago d, [1(1 [H| ; more generally, 
for arbitrary Bravais lattices with general J+j , it was 
shown [T2J that the corresponding spiral, x cos q • n + 
y sin q ■ n, minimizes the classical Heisenberg energy for 
the appropriate wave vector q. See Q for a recent review. 
In the present case, the cluster method provides an al- 
ternate proof (alternative to the Luttinger-Tisza method 
used in [J,[I3]) for the purely Heisenberg case. Because of 
the isotropy of the biquadratic terms, the cluster method 
accomplishes the proof just as easily. 

I list the energies for the various stationary solutions 



hferro = h c (0, 0) = -1 - <Z + 7 

Kudd = K(0, 7r) = -a - 7 
hspiral = h c (9o,—9o) = —j — 
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4(27 - a) ' 

The spiral energy holds only for the condition in 



(10) 




FIG. 2: Phase diagram: a vs. 2 7 

Equating these energies in pairs yields the boundaries 
of the regions shown in FIG. | As a check, to make 
sure no stationary points were missed, I calculated the 
energy difference across boundaries over a mesh of values 
of 9 and 9' varying independently. E.g. I calculated 
h c (9, 9') - Kudd at (27, a) = (1.5, 0.1) and (1.5,0.25), the 
former being a point in the spiral region, the latter in 
the uudd region. The former case showed some negative 
values, the latter only positive values, as must be if the 
phase diagram is correct. 

FIG. [3] shows the variation of q with 7 for a = 0.2. In 
the ferromagnetic and spiral regions, q = 9q, the spiral 
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FIG. 3: q vs 27 at a — 0.2; period or wavelength = 2ir/q 



wave vector; in the uudd region, the significance of q is 
that 2-K/q is the repeat distance of the spin state. If one 
moves from inside the spiral region across its boundaries, 
a lock-in transition occurs at the boundary, seen in the 
special case of FIG. El Interest in this arises because it 
was thought that such lock-in phenomena were caused 
by magnetoelastic couplings or anisotropy, i.e. it results 
from dependence of spins on the lattice [13| • While that 
is probably true in some cases, the present results indi- 
cate another possible cause. Also, as seen in FIG. El the 
transition across the spiral- ferro boundary is continuous, 
while the spiral — > uudd transition is discontinuous. 

In FIG. [51 at a = it is seen that the ferro — > spi- 
ral transition occurs at the well-known value 7 = 1/4. 
At 7 = 0, the transition ferro — > spiral occurs at 
a = —1/2, in agreement with the finding of Thorpe and 
Blume (TB) [8| (their J2 is my —a); in that work only 
a < 0,7 = is considered. Furthermore, they find that 
the state on the line a < —1/2 is disordered; this is not 
inconsistent with the present finding, which implies only 
a spiral in the limit 7 — -> 0; at 7 = the state is indeed 
highly degenerate, since it depends only on the angle be- 
tween nearest neighbors, so that for a given spin S„, S ra +i 
can lie anywhere on a cone with S n as axis and 1/2-angle 
9o, giving a (1-dimensionally) macroscopic entropy. In- 
troduction of the 2nd neighbor Heisenberg interaction 
removes this degeneracy. [14J 

The ferromagnetic transition at a = —1/2 on the line 
7 = shows the following interesting effect: Starting 
from a = 0, adding the extra interaction (the biquadratic 
terms) of sufficient strength causes the transition ferro- 
magnet — > TB disordered state. This is like the inverse 
of the "order-from-disorder" effect 15], which is: in- 
creasing a commonly thought-to-be disordering parame- 
ter, e.g temperature or impurity concentration, can cause 
an entropy reduction. In the present case, adding the bi- 
quadratic terms increases the entropy as a passes through 
- 0.5. I.e. the introduction of an additional interaction 
(usually thought to remove degeneracy, in the spirit of 
the Nernst "theorem"), causes the opposite effect, an in- 
crease in entropy: a "disorder-from-order" effect. 

A surprise is that the spiral state continues for negative 



7. The straight-line ferro-spiral boundary, a — 2j — 1/2, 
continues to — 00 as 7 — » — 00. g or flo vs 7 at fixed 
a < —1/2 changes continuously to zero as the ferro-spiral 
boundary is approached from the right. Nothing special 
happens at 7 = 0, despite the macroscopic degeneracy 
at (and only at) that point. The spiral in this region is 
caused by the competition between the all-ferromagnetic 
Heisenberg exchange and the biquadratic exchange, the 
latter "likes" non-collinear spins with angle between nn. 
spins of 7r/2. The nnn. interaction 7 removes the macro- 
scopic degeneracy (as for the antiferromagnetic case). 

For large positive a and large 7 one sees the up-up- 
down-down phase. This is intuitively reasonable: as al- 
ready discussed, the biquadratic terms with a > are 
very much like the non-interacting Ising model. 

This finding is relevant to the paper, Kimura et al [l6l |. 
which studied the frustrated classical Heisenberg model 
on a square lattice with nearest-neighbor ferromagnetic 
interactions J\ and a 2nd-neighbor antiferromagnetic in- 
teraction J 2 along one diagonal, (1,1) of the square unit 
cell. They were seeking the origin of the "up-up-down- 
down" spin state found in certain manganites, in partic- 
ular HoMnOa. This state shows spin stripes in the a-b 
plane lying along (1,-1), varying up, up, down, down as one 
moves along the (1,1) direction. They presented a phase 
diagram that showed this state at IJ2/J1I > 1. However 
it was recently noted that the correct solution of the as- 
sumed model is quite different, the uudd state occurring 
only in the limit IJ2/./1I —* 00, where it is degenerate 
with a spiral with a 90° turn-angle, propagation vector q 
in the (1,1) direction. [TtJ This realization continued the 
question as to the source of this state, and motivated the 
present study. 

In this connection, one should note another path to the 
uudd state, namely the very different model where the 
nearest neighbor exchange varies from ferromagnetic to 
antiferromagnetic, in continuing periodic fashion. This, 
with no other interactions, trivially leads to the uudd 
state. This one dimensional model is very close to the 
mechanism proposed by Zhou and Goodenough [l8| for 
the same manganites discussed in (l6| . The alternating 
sign of the nearest neighbor exchange interaction in the 
a-b plane of these materials is argued, quite reasonably, 
as being caused by the complex structure of the Jahn- 
Teller distortion. 

I mention two other related works. Girardeau and 
Popovic-Bozic considered the quantum version of 
the model of Thorpe and Blume [|| , showing that in the 
mean field approximation the biquadratic terms are not 
equivalent to replacement by classical spins (unlike the 
Heisenberg terms). Their qualitative conclusions are like 
those of [8|, particularly with respect to the disordered 
phase. Perhaps the earliest paper on the model of com- 
bined Heisenberg-biquadratic interactions on a lattice is 
that due to Rodbell et. al. [2(j for rock-salt structure 
antiferromagnets, MnO , NiO. They assumed negative 
sign for the biquadratic terms (a > in my notation), 
and found large stiffening of the sublattice magnetiza- 
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tion (the stiffening expected for this sign), and within 
their approximation, good agreement with experiment. 
Their model is in a sense closer to the present work in 
that in these structures there are competing Heisenberg 
exchange interactions; however, in these cases the inter- 
actions are consistent with collinear spin states [211 ]. so 
the qualitative behavior is not similar to that found in 
the present work. 

In summary, the ground state of the frustrated Heisen- 
berg model plus biquadratic exchange interactions on a 
linear chain has been solved analytically through an ex- 
act cluster method [||. The case where the Heisenberg 
interactions are all ferromagnetic (unfrustrated Heisen- 
berg model) is included. The phase diagram shows ferro- 
magnetic, spiral, and up-up-down-down spin structures. 
In the unfrustrated Heisenberg case, the spiral is caused 
by the competition between the Heisenberg and the bi- 



quadratic interactions. In this isotropic model, the peri- 
odicity of the spin structures shows a lock-in transforma- 
tion at the boundary of the spiral phase. 

The finite temperature behavior of this model, ex- 
tended to 3 dimensions via the scheme used in the 
ANNNI model 0], (in order to allow long range ordered 
structures), will probably show interesting novel behav- 
ior. As already noted, the ground state in this higher- 
dimensional case is solved by the present results. Also, 
I expect that extension of the ground-state problem to 
higher-dimensional simple cubic models with Heisenberg 
interactions analogous to those in the model of 0, [ItJ , 
and n.n. biquadratic exchange, should again be tractable 
via the exact cluster method [3j . 

I thank S. D. Bhanu Mahanti and Phil Duxbury for 
helpful discussions, and Alex Kamenev and Mark Dyk- 
man for encouragement. 
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